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8 ! ABSTRACT 

r$ Aims. We study the formation of water and methanol in the dense cloud conditions to find the dependence of its production rate on the binding 

1 1 energies, reaction mechanisms, temperatures, and grain site number. We wish to find the effective grain surface area available for chemical 

reaction and the effective recombination timescales as functions of grain and gas parameters. 

^ , Methods. We used a Monte Carlo simulation to follow the chemical processes occurring on the grain surface. We carried out the simulations 
^— ^ ■ on the Olivine grains of different sizes, temperatures, gas phase abundances and different reaction mechanisms. We consider H, O, and CO as 
^j- ' the accreting species from the gas phase and allow ten chemical reactions among them on the grains. 

[-i. ■ Results. We find that the formation rate of various molecules is strongly dependent on the binding energies. When the binding energies are 
■ high, it is very difficult to produce significant amounts of the molecular species. Instead, the grain is found to be full of atomic species. The 
production rates are found to depend on the number density in the gas phase. When the density is high, the production of various molecules 
on the grains is small as grain sites are quickly filled up by atomic species. If both the Eley-Rideal and Langmuir-Hinselwood mechanisms 
qq ' are considered, then the production rates are at this maximum and the grains are filled up relatively faster. Thus, if allowed, the Eley-Rideal 
, mechanism can also play a major role and more so when the grain is full of immobile species. We show that the concept of the effective grain 
• • ■ surface area, which we introduced in our earlier work, plays a significant role in grain chemistry. 
, Conclusions. We compute the abundance of water and methanol and show that the results strongly depend on the density and composition in the 
, gas phase, as well as various grain parameters. In the rate equation, it is generally assumed that the recombination efficiencies are independent 
5_j ■ of the grain parameters, and the surface coverage. Presently, our computed parameter a for each product is found to depend on the accretion 
03 rate, the grain parameters and the surface coverage of the grain. We compare our results obtained from the rate equation and the one from the 
effective rate equation, which includes a. A comparison of our results with the observed abundance shows very good agreement. 
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1. Introduction 

Formation of complex organic molecules in the interstellar 
medium (ISM) is an active subject of research. Out of a host of 
organic and inorganic molecules that have been observed, wa- 
ter and methanol are most certainly two very important organic 
species found both in the gas and solid phases of ISM. The 
abundance of these species in the various regions of ISM is also 
different. It was thought for a long time that the water is one of 
the possible reservoirs of elemental oxygen in the gas phase, 
but recent Submillimeter Wave Astronomy Satellite (SWAS) 
observation found surprisingly a low abundance of water. Snell 
et al. (2000) find that the abundance of water relative to H2 in 
Orion and M17 cloud is between 10~ 10 to 8 x 10~ ll) . The wa- 
ter abundance in the hot cores range from 10~ 6 to 10~ 4 (van 

Send offprint requests to: Prof. Sandip K. Chakrabarti 



Dishoeck & Helmich 1996, Helmich et al. 1996; Boogert & 
Ehrenfreund, 2004). The abundance of water on grains with re- 
spect to the total H column density is typically 10 4 and is the 
most abundant component (Tielens et al. 1991). Similarly, in- 
terstellar methanol has three types of abundance profile: flat 
profiles at CH3OH/H2 ~ 10~ 9 for the coldest sources, pro- 
files with a jump in its abundance from ~ 10~ 9 to ~ 10~ 7 for 
the warmer sources, and flat profiles at a ~ few 10 -8 for the 
hot cores (van der Tak et al. 2000). On the grain surface, the 
methanol abundance varies from 5% to 30% with respect to 
H2O. In some sources, such as SgrA and Elias 16, the abun- 
dance is even less (Gibb at al., 2000). The observed abundance 
for methanol along the line of sight towards high mass and low- 
mass proto-stars is between 0.2 - 2 x 10~ 5 (Gibb et al. 2004; 
Pontoppidan et al. 2003, 2004). Much higher methanol abun- 
dances are found to be associated with the outflows in the re- 
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Model 1 



Model 2 



Model 3 



Fig. 1. Cartoon diagram to show the different reaction schemes 
(C1-C4) that are considered in our calculation. Black and 
shaded circles are representatives of reactive species and the 
clear circles are for nonreactive species. A circle around two 
such species indicates that the new species is created due to 
Eley-Rideal scheme. 



gions of low mass star formation, LI 157-MM and NGC1333- 
IRAS2 2 x 1(T 5 and 2 x 1(T 6 (Bachiller & Perez Gutierrez, 
1997; Bachiller et al. 1998). In other words, either on the grain 
surfaces or in the hotter region, abundances of theses species 
are high. This correlation suggests that these species perhaps 
originate from grains and their productions in the gas phase are 
inadequate. Therefore, the understanding of the formation of 
water and methanol on grain surfaces is of primary importance. 

The grain surface reactions were first introduced to explain 
the formation of molecular hydrogen (Hollenbach & Salpeter, 
1970). Since then it has been used very extensively by sev- 
eral authors (Watson & Salpeter 1972ab; Allen & Robinson 
1975, 1976, 1977; Tielens & Hagen 1982; Hasegawa & Herbst 
1992; Charnley 2001; Stantcheva et al. 2002, Green et al. 
2001; Biham et al 2001; Stantcheva et al 2002). These stud- 
ies mainly belong to two different categories, the determinis- 
tic approach and the stochastic approach. In the determinis- 
tic approach, one can completely determine the time evolu- 
tion of the system once the initial conditions are known. The 
rate equation method belongs to this category. This method is 
very extensively used by several authors to study the grain sur- 
face chemistry (Hasegawa & Herbst 1992; Roberts et al 2002; 
Acharyya et al. 2005). However, this method is only applicable 
when there are large numbers of reactants on the grain surface. 
Given that the interstellar medium is very dilute, very often 
this criteria is not fulfilled and this method cannot be applied. 
But this method is computationally faster and can very easily 
be coupled with the gas phase reactions. In the stochastic ap- 
proach, fluctuations in the surface abundance due to the statis- 
tical nature of the grain is preserved. The Monte Carlo method 
and the Master equation methods belong to this category. Both 
these methods are used by several authors (Charnley 2001; 
Stantcheva et al 2002, Green et al. 2001; Biham et al 2001; 
Stantcheva et al. 2002 ). Its major disadvantage is that it takes 
enormous computational time. Coupling of the Monte Carlo 
method (for grain surface reactions) and rate equation method 
(for gas phase reactions) is extremely difficult. Although the 
Master equation method can be coupled with the rate equa- 
tions, but it is disadvantageous because one has to solve a large 
number (ideally infinite) of reactions. 
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Fig. 2. Evolutions of the number of a few selected species on 
a grain having 10 4 sites are shown with respect to the surface 
coverage of a mono-layer (along X-axis), for various sets of 
binding energies (Models 1-3) and for different mechanisms 
(C1-C4). These simulations were carried out for high abun- 
dances of the accreting species (Table 1). Different styles of the 
curves are marked with the species names on the right. Note 
that methanol and CO2 are absent in the Model 1 simulation 
where the binding energy was very high. 



Recently, Chang et al. (2005) have argued that the stochas- 
tic methods used so far can also lead to error because the rate 
of reaction is determined by the rate of hopping (or tunnel- 
ing) of a hydrogen atom from one site to the nearest neigh- 
boring site multiplied by the probability of finding a reactant 
partner in this site. This is also an average treatment since, 
on any given grain, the reactant partner is unlikely to lie in 
the nearest-neighboring site. They used a continuous random 
work technique to study the formation of molecular hydro- 
gen. Chakrabarti et al. (2006ab) used a similar method that 
keeps track of each individual reactant and their movements 
and calculated the effective grain surface area involved in the 
formation of molecular hydrogen in the interstellar clouds. 
Chakrabarti et al. (2006, hereafter referred as Paper-I) illus- 
trated that the formation rate per unit grain site itself strongly 
depends on the nature and size (i.e., no. of sites) of the grains 
for a given gas phase condition (abundance, temperature, etc.). 
In Paper-I, this was demonstrated only for H forming H2 
molecules. In the present paper, we carry out a similar anal- 
ysis where we consider the accretion of H, O, and CO onto the 
grain surface and show that the formation rates of water and 
methanol are indeed dependent on the intrinsic and extrinsic 
parameters of the system. The plan of this paper is the follow- 
ing. In the next section, we discuss various mechanisms of re- 
actions on the surface. In Section 3, we discuss the procedure 
of our computation. In Section 4, we describe various models 
used for our simulations. In Section 5, we present our results. 
Finally, in Section 6, we draw our conclusions. 
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Table 1. Gas-phase abundances used 
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Fig. 3. Same as in Fig. 2, but the low abundances of the accret- 
ing species have been chosen. 



2. Mechanisms of reactions on grain surfaces 

A thorough understanding of the surface reaction mechanisms 
requires knowing of the basic physical processes involved 
when the gas phase atoms and molecules interact with the 
grains. The first step is accretion, i.e., landing various species 
onto a grain surface. In our case, only H, O, and CO are taken 
as the accreting species onto the grain surface. In the next step, 
the accreted species will react to form various new species. 
There are two reaction schemes, the Langmuir-Hinselwood 
(LH) mechanism and the Eley-Rideal (ER) mechanism. In the 
LH scheme, the gas phase species accretes onto a grain and 
becomes equilibrated with the surface before it reacts with an- 
other atom or molecule, and in the ER reaction scheme, the 
incident gas phase species collides directly with an adsorbed 
species on the surface and reacts with that species. In such a 
mechanism, the reactant generally does not become trapped at 
the surface and it is unlikely to be sensitive to the surface tem- 
perature (Farebrother et al. 2000). In our study, we considered 
that the reaction on a surface can occur through both LH and 
ER mechanisms. However, we assume that the molecules are 
trapped after reaction due to their high binding energy. 

The binding energy of the incoming species strongly de- 
pends on the species itself and on the way the interaction pro- 
ceeds. The incoming species might get trapped in a shallow 
potential well in a physisorbed site. The interaction is mainly 
due to mutually induced dipole moments or it might also form 
a strong covalent bond. Recent studies have found evidence of 
both physisorption and chemisorption processes taking place 
on a grain surface. However, for the chemisorption, a higher 
kinetic energy is involved, so this type of interaction is not rel- 
evant except in very special astrophysical conditions. We have 
considered only weakly bound species, i.e., physisorbed atoms 
and molecules. The typical energy for physisorption is around 
0.1 eV or ~ 1000 K. If Ed denotes the binding energy for 
physical adsorption and E/, the potential energy barrier, then 
Ed must be greater than Ej, for the species to diffuse from one 
site to the other. 



Species 


high (cm 3 ) 


low (cm i ) 


H 


1.10 


1.15 


O 


7.0 


0.09 


CO 


7.5 


0.075 



For a chemical reaction to occur, an accreted species has 
to scan the grain surface in search of a reaction partner. There 
are two physical processes that can provide the mobility for the 
accreted gas phase species, the thermal hopping and tunneling. 
As mentioned earlier, Hollenbach & Salpeter first introduced 
the grain surface chemistry to explain the high abundance of 
molecular hydrogen. They assumed that, within the grains, the 
hydrogen atoms move from site to site by a quantum mechani- 
cal tunneling process. But from the recent experimental results 
of molecular hydrogen formation, it was found that the mo- 
bility of hydrogen on grains is primarily due to the thermal 
hopping (Pirrenello 1997ab, 1999). The tunneling time very 
dependent on the mass of the particle and the barrier thick- 
ness, a unknown parameter. Therefore, it is not widely used in 
the astrophysical models. The classical papers like Tielens & 
Hagen (1982) and Hasegawa & Herbst (1992) considered the 
tunneling for hydrogen and thermal hopping for other simple 
atoms and molecules. We considered both tunneling and ther- 
mal hopping for the hydrogen atoms and only thermal hopping 
for others. 

If we consider the time scale for hopping and tunneling for 
the hydrogen atom and hopping for other species, we observe 
that the mobility of hydrogen atom is much greater than for 
the other reactive species. The hopping and tunneling time for 
hydrogen is ~ 7.4 x 10~ 09 and 2.0 x lO^'s, and the hopping 
times for atomic species like O and CO are ~ 2.4 x 10~ 2 s and 
~ 5589s , respectively. Therefore, the hydrogenation reactions 
are the dominant reactions on grains. 

2.1. Accretion rate 

Accretion is the process by which grain receives matter. We 
define the accretion rate [r acc (i)] for a given neutral species i in 
the units of s _I , as 



c (0 = Si<TVin gas (i), 



(1) 



where s,- is the sticking coefficient (taken as 1 for all the three 
species), v, the velocity (cm s _I ), n gas (i) is the gas phase num- 
ber density (cm 4 ) of the z'-th species, and cr the grain cross- 
section (cm 2 ). In this paper, we carry out our calculation usu- 
ally for two different sets of gas phase number densities, ex- 
cepts towards the end where we compared results from sev- 
eral number densities. Thus, we have two sets of accretion 
rates. Number densities are taken from Stantcheva et al. (2002). 
These are listed in Table Q] 

2.2. Diffusion rate coefficients and binding energies 

The diffusion rate coefficients (k) are the sums of the rates fi, 
(s -1 ) the reactive species need to traverse an entire grain (tdiff 
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Fig. 4. Rate of production of H2O is shown against the sur- 
face coverage for the four different schemes mentioned in the 
text. The production is high when the Eley-Rideal mechanism 
is taken into account, while the production is lower when the 
reaction takes place at the next time step. The rate is computed 
using (a) the total production and total time taken at a given 
instant, (b) the change in the number of H2O for a fixed change 
in parcentage of the surface coverage. 



is the diffusion time). In the present paper, we are interested 
in showing the size dependence of the recombination process. 
Thus, we choose a grain to contain 10 4 to 10 5 sites. The rate 
coefficient is usually multiplied by a factor k that accounts for 
any non zero activation barrier (Hasegawa et al. 1992). 

The rates depend strongly on E/,, the barrier energy. We 
considered three different sets of barrier energies. In the first 
set, we used the binding energies that come from the ear- 
lier works (Allen & Robinson, 1977; Tielens & Allamandola, 
1987; Hasegawa & Herbst, 1993). We show them in Table 3. 
In this set, we considered that the hydrogen diffusion is caused 
by thermal hopping. Thus we ignore tunneling. In the second 
set, we used the same set of binding energies except for hy- 
drogen diffusion procedure. Instead of hydrogen diffusion, we 
considered tunneling. Since the hydrogen diffuses much faster 
on grains than any other species, hydrogenation reaction is the 
dominant reaction on grains. This is why we did done our simu- 
lation for both thermal hopping and tunneling of hydrogen. The 
third set is based on recent findings of Pirronello et al. (1997, 
1999) as interpreted by Katz et al. (1999), which show that 
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Fig. 5. Rate of production of a few selected species as a func- 
tion of surface coverage at different temperatures shown for a 
grain having 10 4 sites for Model 2 and C2 scheme. These sim- 
ulations were carried out for low abundances of the accreting 
species. 
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Fig. 6. Plot of the time to build a mono-layer as a function of 
temperature for a grain having 10 4 sites and for Model 2 by 
considering the C2 scheme. These simulations were carried out 
for low abundances of the accreting species. 

atomic hydrogen moves much more slowly than what is used 
in various simulations. Therefore, in this set we used hydro- 
gen hopping rates from Katz et al. (1999) and barrier energies 
for other species are increased proportionately (Table 2). This 
means that the ratio between the barrier energies of these two 
tables for different species are exactly the same as for H (i.e., 
2.87 for Ei, and 1 .066 for Ej). The ratio between the desorption 
energies of Table 2 and 3 was also kept fixed for simplicity. We 
chose E/, = O.3E4 in Table 3 for all species except for H. 

Using these three sets we have constructed three models. 
The results are shown according to the increasing diffusion rate 
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Table 2. Energy barriers against diffusion and desorption in K 
for the olivine grain (from experiment) 



opcuics 


F, 


F , 


T-[ 


287 


373 


o 


689 


853 


OH 


1085 


1343 


H 2 


387 


479 


o 2 


1042 


1290 


H 2 


1601 


1982 


CO 


1042 


1290 


HCO 


1300 


1609 


H 2 CO 


1515 


1876 


CH3O 


1868 


2313 


CH3OH 


1774 


2195 


co 2 


2153 


2664 



of hydrogen. Thus in Model 1, we used the slowest diffusion 
rate of hydrogen and hence the binding energy of Table 2. In 
Model 2, we used the binding energy from Table 3 (below). 
Finally in Model 3, the hydrogen tunneling is considered and 
the binding energy is taken from Table 3 for other species. 

Table 3. Energy barriers against diffusion and desorption in K 
for the olivine grain (See text for references.) 



Species 


E b 


E d 


H 


100 


350 


O 


240 


800 


OH 


378 


1260 


H 2 


135 


450 


o 2 


363 


1210 


H 2 


558 


1860 


CO 


363 


1210 


HCO 


453 


1510 


H 2 CO 


528 


1760 


CH3O 


651 


2170 


CH3OH 


618 


2060 


co 2 


750 


2500 



3. Simulation procedure 

We used a Monte-Carlo simulation method of studing the 
chemical evolution on grains in the presence of H, O, and CO 
as accreting species. For the sake of simplicity we assume that a 
grain is square shaped having S number of sites. For concrete- 
ness we chose olivine grains having the surface density of sites 
s = 1.5 x 10 15 crrT 2 (Hasegawa et al. 1992). We assumed fur- 
ther that each site has four nearest neighbors, as in an fcc[100] 
plane. To mimic the spherical grain structure, we assumed a pe- 
riodic boundary condition, i.e., a species leaving the boundary 
site of a grain on one side enters back into the grain from the 
opposite side. 

In our method of simulation, self-crossing paths of ran- 
domly walking atom are automatically included (some times 
called the back-diffusion in the literature, e.g., Chang et al. 
2005). The minimum time step of our simulation is assumed 
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Fig. 7. Variations in jH 2 o(t) with surface coverage for different 
reaction schemes for Model 2. Note that when the Eley-Rideal 
scheme is included (CI and C4) the production efficiency in- 
creases rapidly as the grain is filled up. 

to be the time taken for hopping or tunneling (whichever is 
smaller) of a H atom from site to site. If r acc is the astro- 
physically relevant accretion rate under consideration, the z'-th 
species will have to be dropped after every l/r acc seconds. The 
location of the accreting ;-th species is obtained from a pair 
of random numbers (R x , R y ?>, themselves obtained by a random 
number generator. This pair would place the incoming species 
at (J, k)-th site of the grain, where j and k are the nearest in- 
tegers obtained using Int function: j = int(R x * n + 0.5) and 
k = int(R y *n+0.5). Here, n = y[S, being the number of sites on 
a square. Now during the hopping or tunneling process, when 
one species enters the site that is already occupied by another 
reactant surface species, it will form a molecule if activation 
barrier energy permits. Thermal evaporation of a species is also 
handled by random numbers. After each hopping and tunneling 
time, we generate a random number (R t ) for ;'-th species, if R, is 
less than WJah, or Wi/a, r Then we allow that species to evap- 
orate. Here, W, is the desorption rate and and a tj are the 
rates of hopping and tunneling of the ;'-th species, respectively. 
We carry out our simulation up to the time within which one 
monolayer is produced. We wish to make a remark in passing 
that a totally different random number generator was also found 
to yield a similar result, so thus we believe that the results we 
present in this paper do not depend on a random number gen- 
erator. 

Keeping in mind that there are two types of reaction 
schemes, the landing of a gas phase atom or molecule can have 
one of these two distinct fates. It can land either on a vacant site 
or on an occupied site. When the landing is on a vacant site, 
the follow-up computation is straightforward. The species will 
simply start scanning the surface for a reactive partner. When 
it lands on an occupied site, however, there could be two possi- 
bilities, (i) The incoming species can land directly on a reactive 
species but it may or may not combine with it. If not, it will look 
for an alternate site for itself, (ii) The incoming species can land 
directly on a nonreactive species and look for an alternate site. 
However, looking for an alternate site itself could be achieved 
in two ways: in the same time step (which is more likely) or 
in the next time step. Thus we can have a total of four differ- 
ent models in the simulation. These are illustrated in a cartoon 
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Fig. 8. Variation in y^oit) with the surface coverage for differ- 
ent grain sizes. Here, the binding energies of Model 2 are used. 
Note that jH^o(f) is decreased with increasing grain size since 
the probability of capturing an incoming species increases with 
the grain size. Here C2 scheme has been chosen. 
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Fig. 9. Variation in a# 2 o(f) with surface coverage for different 
reaction schemes as in Fig. 7. The formation rate of the new 
species goes down (i.e., aH 2 o(t) goes up) for C2 and C3 since 
the species is created in the next time step. 



diagram C1-C4 in Fig.Q] In CI, both the LH and ER schemes 
are considered, i.e., the incoming atom/molecule is allowed to 
land on both the unoccupied and occupied sites. With the reac- 
tive species, it combines to form a new species (marked with a 
common circle). It skips the nonreactive species and looks for a 
new site in the same time step. C4 is otherwise the same as CI, 
but the incoming species looks for a new site at the next time 
step after it lands on a nonreactive species. In C2, the incoming 
species lands on a reactive species, but does not combine and 
insteads finds a new site at the same time step (LH mechanism). 
C3 is the same as C2, but we allow the landing at a new time 
step. 

Table 4. Surface Reactions in the H, O, and CO model 



Number 


Reactions 


E a (K) 


1 


H+H -> H 2 




2 


H+O OH 




3 


H+OH H 2 




4 


H+CO HCO 


2000 


5 


H+HCO -> H 2 CO 




6 


H+H 2 CO -> H3CO 


2000 


7 


H+H 3 CO -» CH 3 OH 




8 


O+O 2 




9 


O+CO -> C0 2 


1000 


10 


O+HCO -> CO,+H 





4. Model simulations 

We consider that only three species, namely, H, O, and CO are 
accreting in the grain phase, and we assume that ten types of 
reactions can go on among these constituents. These are listed 
in Table [4] Let us consider a reaction of the type, 

a + b — > c. (2) 



The rate of variation of any molecular species c with time is 
usually taken as (see, Paper-I and references therein) 

dn c 

— = r ab n a n b /S (3) 
at 

where, n a , n b , and n c are the concentrations of the species a, 
Z?,and c on a grain. Here, r ab is the rate of formation of c out 
of a and b, and S is the number of sites on the grain. The rate 
r ab /S is approximately the inverse of the time required for an 
atom to visit nearly all the adsorption sites on the grain surface. 
This is because, in two dimensions, the number of distinct sites 
visited by a random walker is linearly proportional to the num- 
ber of steps, up to a logarithmic correction. But in a realistic 
situation, this is not the case for two reasons (Paper-I). The 
first is that the number of hoppings to find a reactant partner 
depends on the number of atoms or molecules already occupy- 
ing the grain. When the grain is almost empty, a species has to 
hop several times, and when the grain is full, it requires less 
of hops, so the 'size' of a grain is relative as far as the hop- 
ing element is concerned. The second reason is the blocking 
effect. This phenomenon is seen very often on a grain having 
a mono-layer. On a grain surface, the interaction of species is 
generally considered to be with those located in the four nearest 
neighboring sites. A possible scenario would be that a reactive 
species is blocked in all the four directions due to the presence 
of various species with which no reaction is allowed. For in- 
stance, if a CO molecule lands on a site surrounded by H2O 
molecules, it cannot react and will wait there itself for its des- 
orption unless the blocked species themselves desorb and free 
the sites for the landed species to move around. The effect be- 
comes more severe as the grains start getting filled up. In the 
latter case, the blocking effect will weaken when the blocked 
element overcomes the energy barrier. As the surface is popu- 
lated, these types of effects become more and more prominent. 
To understand the influence of these processes on the recom- 
bination, we define two parameters, namely, y and a as given 
below. 
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Fig. 10. Same as in Fig. 8 except that a^oit) has been plotted. 
As the grain size gets smaller, a^o increases, i.e., the recom- 
bination efficiency decreases. See text for details. 

If t c is the "average time" required to form a species c af- 
ter the species a and b are simultaneously present in our sim- 
ulation, then -|r ~ n a rib/t c , where t c = r a b/S. The rate with 
which the species c is formed is l/t c . Because of the two ef- 
fects mentioned above, the actual rate may be assumed to be 
r a blS y c instead of r a b/S . From the above equation, 

log[r ab t c ] 

lc = . c ■ (4) 
logS 

As an example, let us consider the calculation of y c for H2O. 
Since H has a very small desorption time scale compared to the 
other species, it will desorb from the grain if it does not find 
a reaction partner within its desorption time scale. If H gets 
another OH to form H2O, then we calculate dtn 2 o by noting 
down the time required to form an H2O after the accretion of 
an H. At the beginning of the simulation, most of the sites are 
empty and one H will have to more or less scan the whole of 
the grain to find its reactant partner, and hence will take longer 
time to form a molecule. On the contrary, when most of the 
sites are filled with the different species, it will take less time 
to form a molecule and y c should decrease. When a molecule 
is formed due to direct hitting, i.e., through the ER mechanism, 
y c is clearly zero, since production rate is independent of the 
population of the species on the grain. Thus, if ER mechanism 
is included, y c tends to become negligible as the grain is popu- 
lated. 

While y is derived from the time interval between the re- 
lease of the reactants and the actual reaction to take place, an- 
other quantity, namely a may be defined to consider the cumu- 
lative effects and thus depends on the total time taken since the 
start of the accretion onto the grain. Because of the two compet- 
ing processes mentioned at the beginning of this section, there 
will be a net change in the rate equation itself with time, since 
a itself will be time dependent. To quantify this, we replace 
r a b/S in Eq.[3]by r a b/S a ^K Thus, the Eq.[3] which is rewritten 
as 

log^ = log[r ab n a n b IS aM l (5) 



can be used to obtain a c (t) as 

a c (t) = log[ ( -!^^-yiog(S). (6) 

dt 

Note that, for reactions between the similar species, their will 
have to be a factor of 1/2 inside the logarithmic expression 
(see, Paper-I for H2 formation) to reduce the possibility of dou- 
ble counting. As the surface is populated by a number of reac- 
tant species, the recombination time scale would vary with time 
and thus a c (f) will also vary. As mentioned before, the time t 
used in computing a is the total time taken since the grain starts 
populating. On the other hand, y is an average quantity that is 
computed instantaneously. 

We run the simulations to show that a{t) and y are not 
necessarily unity. This means that the surface chemistry is not 
just a matter of random hopping, since the degree of random- 
ness is seriously modified with the population and the abun- 
dance of species on the grain surface. The Monte-Carlo simu- 
lation reflects a more realistic computation, but it is very time- 
consuming. In a large-scale computation, it is impractical to do 
such a simulation. Indeed, we also show that only when a as 
given above is used in the effective rate equation (where every 
S is replaced by S a , where a, is different for different species), 
then the abundances of each species obtained from such an 
equation agree with the results of the Monte-Carlo simulation. 
Thus, one may as well solve the effective rate equation using 
a obtained from our calculations. We have made a table for 
a (Table 7 , below) for future. However, discussions of its ap- 
proximate value, shape, and dependence is beyond the scope of 
this paper and will be discussed elsewhere. 

In what follows, we ran twelve different cases. For each set 
of binding energies, all the four CI, C2, C3, and C4 cases are 
considered. Due to the complexity of the system, we restricted 
ourselves only up to the formation of a mono-layer on the grain 
surface. 

5. Results 

We first carried out the simulation for a grain having 10 4 sites. 
We considered two different sets of gas-phase abundances in 
our study, one with a high density and the other with low den- 
sity. For each set of gas-phase abundances, we consider three 
sets of binding energies and for each set of binding energies we 
have four cases, CI, C2, C3, and C4 as described in Section 3. 
Thus for each set of gas-phase abundance, we essentially have 
twelve different cases. 

5.1. Grain surface abundances 

In Fig. |2] we show the results for the high gas-phase abun- 
dances. In the first column, the plots correspond to the binding 
energies of Model 1 as mentioned in the Table|2]are shown. The 
binding energies are high and consequently the diffusion rates 
are low. Let us first consider CI and C2 cases. We note that, re- 
spectively, only 16% and 1 1% of the grain deposition is water. 
This low production of water is due to the slow diffusion rate 
and due to the absence of sufficient number of hydrogen atoms 
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Fig. 11. Variation in a with the surface coverage for various 
molecular species that form on the grain surface. All of them 
deviate from unity. This case is for Model 2 with the C2 
scheme. 
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Fig. 12. Variations in a for water and methanol (averaged over 
the surface coverage) as functions of the grain radius and no. 
of sites on the grain. The filled circles are for H2O and the 
filled boxes are for methanol. The corresponding fitted curves 
are also shown. They are extrapolated to larger grain sizes. 



to react with the oxygen atom. In CI, the excess of produc- 
tion of water by four per cent is mainly due to ER mechanism. 
There is hardly any methanol produced due to slow diffusion 
rate and due to presence of an activation barrier energy. In CI, 
around 19% of molecular oxygen is produced, but its absence 
in C2 suggests that it is produced due to the ER mechanism. 
Around 50% of the grain surface abundance is CO. The time 
required to grow one mono layer for CI and C2 is 1600 years 
and 1300 years, respectively. The grain is very quickly filled 
up by the accreted species. The results for the cases C3 and C4 
differ only slightly in comparison to C2 and CI, respectively. 
Only a small increase in the molecular species and decrease in 
the atomic species is observed. The times required to grow a 
mono layer in C3 and C4 are around 13000 years and 36000 
years, respectively. This longer timescale is due to the fact that 
in these cases atoms/molecules are not bound to accrete on the 



c 



Hopping 




Surface Coverage(%) 

Fig. 13. Variation in the fractional productions of H2O through 
ER and LH schemes against the surface coverage of a mono- 
layer. For low surface coverage, H2O production is dominated 
by the hopping process and when surface coverage is higher, 
the Eley-Rideal mechanism contributes in a major way to- 
wards the H2O production. The uppermost and the lowermost 
curves are for the Model 2 energy and CI scheme, with the low 
abundances of the accreting species. We compare the same for 
Model 1 energy with CI method, with the high abundances for 
the accreting species (middle two curves). 



same time step. In the case of C4 we continue our simulation 
up to the time when all O and OH are converted to H2O, hence 
the time taken is much longer. 

The second and third columns show the plots correspond- 
ing to binding energies as mentioned in Table [3] In the second 
column, the hydrogen diffusion is due to the thermal hopping, 
whereas in the third column this is due to tunneling. These 
binding energies are low and the diffusion rates are conse- 
quently high. In all the four cases, the most dominant species 
are CO, O2, and water. Water abundance is very similar, vary- 
ing between 11% and 19%. A very small amount of methanol 
has been produced. In all the cases, the production of a signif- 
icant amount of molecular oxygen suggests that the ER mech- 
anism may have an active role to play. The time required to 
grow a mono-layer is around 1500 years. The third column also 
shows a similar trend of the results. 

In Fig. [3] we show the results for low gas-phase abundance. 
A significant amount of water is produced in all the cases, 
since the H/O and H/CO is relatively high. The presence of 
the atomic oxygen and carbon monoxide is also much lower, 
especially for Models 2 and 3. This is because, in the last cases 
the mobility of H being many orders of magnitude higher hy- 
drogenation process is more efficient and that uses up O and 
CO to create water and methanol. This makes the production 
of methanol significant for Models 2 and 3 and insignificant 
for Model 1 . Compared to the high abundance cases (Fig. 2), 
the accretion rate for low abundance is much lower, so the grain 
is not has full of atomic species as before. The time required to 
grow a mono-layer is around 0. lMyr for both the cases. A sim- 
ilar trend is found for C3 and C4, but it takes around lMyr to 
fill up the grain. The second and third columns of Fig. 3 show a 
similar result. In all the cases, the most dominant species is wa- 
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Fig. 14. Comparisons of the number of H2O and CH3OH on a 
grain of 10 4 sites as obtained by the rate equation method, the 
Monte-Carlo method, and the effective rate equation method. 
Here we chose the Model 2 energies and C2 with the low abun- 
dances case. The agreement between the Monte-Carlo result 
and the effective rate equation result is very good. 

ter. The abundance of water varies between 46% and 61%. The 
next dominant species is methanol with an abundance between 
26% and 33%. But only a trace amount of molecular oxygen 
is produced. Table [5] shows the percentage abundances for se- 
lected species, and Table [6] shows the time required to grow a 
mono-layer. 

Figure 4(a-b) shows the rate Rh^o of production of H2O 
for various mechanisms. This has been computed in two 
ways. Figure 4a shows the cumulative average rate where 
Rh^o =ihO(t)/t. We clearly see that, unless the production rate 
is very high (e.g. in case of CI), the cumulative rate generally 
falls off very rapidly with time. In the CI scheme, the abun- 
dance of H2O is higher because the formation of H2O is domi- 
nated by direct hitting rather than hopping. In Fig. 4b, we plot 
Rh 2 o - ~JT"> where dH 2 is the production of H 2 in dt time 
in which a fixed percentage dp% of the surface coverage is in- 
creased. Thus, in effect, it is the instantaneous rate. For simplic- 
ity, we choose dp = 10% so as to get smooth curves. In the case 
of the LH scheme, i.e., for C2 and C3, as the surface starts to 
populate, the rejection of the incoming species become impor- 
tant and thus the rate of production decreases with the surface 
coverage. 

5.2. Temperature dependence 

Figure 5 shows the production of various species as a function 
of surface coverage at 10 K, 15 K and 20 K. We carried out 
this computation only for Model 2 and taking low abundances 
of the accreting species. In all the three temperatures, the wa- 
ter is the most abundant species. But CH3OH is only efficiently 
formed in 10 K and 15 K. At 20 K, methanol production is 
significantly low because the desorption time scale of atomic 
hydrogen is much shorter and, at the same time, the activation 
barrier energy is much higher to produce methanol. It is in- 
structive to study the time taken to fill up the grain. In Fig. 6 
we plot the time required to build a mono-layer on a grain of 



10 4 sites as a function of the temperature of the grain. At higher 
temperatures, the time taken is lower because the grain is full 
of nonreactive CO (Fig. 6). 

5.3. Variation in y and a 

In Fig. 7, the variation in 'catalytic capacity' for water, namely, 
7H 2 o(t) (Eq. 4) is shown as a function of the surface cover- 
age. The figure shows that 7h 2 o(0 gradually decreases with 
the surface coverage. This is because an increase in the sur- 
face coverage will decrease the number of vacant sites and will 
decrease the time for searching its reactant partner to recom- 
bine. The blocking effect also shows its effect by reducing jh 2 o 
very rapidly. The lower the value of y#,o, the more efficient the 
recombination. This parameter is dependent on the number of 
sites on the grain. Figure 8 shows the change of jh 2 o with the 
number of grain sites. Naively, one would have expected that 
jHiO to remain unity always. In that case, jh 2 o would be in- 
dependent of number of sites in the grain. However, we find 
that jh 2 o goes down with increase in S . From Eq. (4) we note 
that, when S is increased from 10 4 to 10 5 , the value of jh 2 o 
should go down by about 20%. This is roughly what is seen in 
the figure. 

In Fig. 9, the variation in aH 2 o (Eq. 6) for water is shown 
as a function of time, which is represented here as the surface 
coverage. Here 10 4 sites are chosen. This exponent deviates 
significantly from 1. This deviation is mainly due to blocking 
the reacting species by the nonreactive species and the popula- 
tion of various species on the grain surface. This is particularly 
true for cases C2 and C3, since the ER mechanism was not 
assumed. For those with an ER mechanism (CI and C4), the 
numbers actually go down, so the formation rate goes up. In 
Fig. 10, we plot the variation in an 2 o as a function of surface 
coverage for a different number of sites on the grain. As the 
number of sites increases, a^o is seen to decrease as in the 
case of jh^o- We thus observe that the production rate of any 
species goes up for larger grains. To judge the implication of 
this in a cloud, we note that the dependence of the number den- 
sity on the radius r goes typically as r~ 35 (Mathis et al. 1977), 
while the dependence of the production rate is r" with n ~ 1 . 
Thus the production is still dominated by the smaller grains. 
However, smaller grains are lesser stable and could easily evap- 
orate. This aspect clearly requires more thorough study. 

In Fig. 1 1, we plot a for various species as a function of sur- 
face coverage when 10 4 sites and the Model 2 with C2 mech- 
anism were chosen. The most important observation from this 
is that a is not unity for any of these species. Second, up to 
about 30% of the surface coverage a roughly remains constant, 
but after that all of them are going up monotonically as the sur- 
face gets filled up. Here the ER scheme is not used, and the 
formation rate goes down with coverage. 

Since we generally see that a deviates from unity for 
smaller grains, it is pertinent to ask at, roughly what grain sizes 
a becomes unity. This would give us a region in which the clas- 
sical rate equation would be valid. In Fig. 12, we show the vari- 
ations in a for water and methanol are shown as functions of the 
grain radius and number of sites on the grain. The correspond- 
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Fig. 15. The ratio of methanol and water abundances on the 
grain surface plotted against the number density of the accret- 
ing gas. From the figure it is clear that as the number density 
goes up, methanol production goes down. The observed abun- 
dance of solid state methanol w.r.t. water is 5%-30%. From the 
figure it is to be noted that the region lies in the range 2.6 x 10 4 
to 8 x 10 4 . These simulation were carried out for Model 2 by 
considering the C2 method. 



ing fitted curves are also shown. They are extrapolated to larger 
grain sizes. We note that, a ~ 1 at around r ~ 3 - 5yum. Since 
the computational time to actually carry out the simulation at 
this size is prohibitively high, it could not be independently ver- 
ified if this conclusion is rigorously valid. On the other hand, 
the number of grains at this size is expected to be negligible, 
and thus the exact knowledge of this limit may not be essential. 

5.4. A comparison between LH and ER schemes 

In Fig. 13, we plot the fractional production of H2O for two 
different schemes at a given surface coverage to show the rel- 
ative importance between the ER and LH schemes. Initially, 
when most of the sites are empty, the production of H2O is 
mainly due to the hopping mechanism i.e., the LH scheme. 
As the grain starts to populate, the ER scheme starts to play 
an important role. In the figure, the uppermost and the lower- 
most curves are obtained using the low abundance set with CI 
model. The other pair of curves is for a higher abundance set 
in the gas phase. Here the effect is more prominent. In general, 
the abundance of all the species is always a few % higher if the 
ER scheme is allowed. 

5.5. Comparison of results with the effective rate 
equations 

It may be instructive to check if the rate equations are modified 
with IS replaced by r a t,/S a everywhere, and they are solved 
by the usual method (Acharyya et al. 2005), then whether the 
results become comparable to those obtained from our Monte- 
Carlo method. In Fig. 14, we made a comparison of the abun- 
dances of two of the important species, namely, H2O and 
CH3OH in different methods. The results of our simulation 



clearly agree with those from the effective rate equation very 
well, but not with the original rate equation. This indicates that 
the as we introduced are very important and their effects are to 
be taken into account for accurately estimating of the species 
from the grain chemistry. Indeed, using the well-known MRN 
model of grain size, we find that water and methanol will be 
underproduced by ~ 25% and ~ 33%, respectively, when the 
effective rate equations are used instead of the usual rate equa- 
tion. 

Since a depends on both the gas and the grain parameters 
in a non-trivial way, it is difficult to obtain any analytical for- 
mula or fitted curves without repeating the simulation many 
times. This is beyond the scope of this paper and will be re- 
ported in future. Nevertheless, we provide in table 7, a list of a 
for olivine grains at 10K for Model 2 binding energies and C2 
scheme of interaction as a function of the surface coverage and 
number of sites. For other sites and coverages, one can use an 
interpolation method. We do not give any table for y since it is 
not immediately used in the rate equation. 

6. Concluding remarks 

In this paper, we studied the formation of water, methanol, and 
other related species on a grain surface using a Monte-Carlo 
method. We used three different sets of binding energies, two 
different gas-phase abundances and considered both the ER and 
LH mechanisms. We considered all the four possible ways (Cl- 
C4) of the landing of a species on the occupied site of a grain. 
Besides the formation of these species, we also calculated two 
parameters y and a that represent the recombination efficien- 
cies to show that they are indeed dependent on the number of 
sites on the grain and the populations of various species on a 
grain surface. 

We found that the formation of various molecules is de- 
pendent on the binding energies. We find that, when the higher 
binding energies are used, it is very difficult to produce a sig- 
nificant amount of the molecular species, instead, the grain 
is found to be full of atomic species. The formation of these 
species is also dependent on the gas phase density. We found 
that, for the high density case, the production of various 
molecules is small and the grains are filled up relatively quickly 
by atomic species. We found that, if both the ER and LH mech- 
anisms are considered, then the production is always high and 
the grain is filled up very quickly. As expected, we find that, 
when the grain is more or less empty, the LH scheme is most 
important. The ER scheme starts dominating more and more as 
the grains are filled up. We have verified that our results from 
the Monte-Carlo simulation match with that from the effective 
rate equation. 

It is important to see if we can put a constraint on the model 
parameters from the observational results. We already men- 
tioned in the Introduction that the abundance of water relative 
to H2 in cold clouds is between 10~ 10 and 8 x 10~ 10 and be- 
tween 10~ 6 to 10~ 4 in hot cores. The abundance of methanol 
with respect to H2 is ~ 10~ 9 for coldest clouds and between 
~ 10~ 9 to 10" 7 for warmer clouds and ~ a few xl0~ 8 in hot 
cores. In the grain surfaces, the solid state methanol abundance 
should be at the most 30% with respect to H2O. If we glance at 
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Table 5 and assume that the surface coverage on multiple lay- 
ered grains is similar to that on a mono-layer, then it is clear 
that the case of Model 1 and high abundance cases for Models 
2 and 3 are not very relevant for the production of methanol. 
Indeed, the low abundance cases with Models 2 and 3 produce 
a similar solid state methanol with respect to H2O (~ 43% - 
~ 57%). Corresponding high abundance models produce only 
4% -6%. Thus the observed abundance of ~ 30% must be from 
some intermediate region of the cloud. In Fig. 15 we show the 
ratio of methanol and water as a function of the number den- 
sity of the gas phase, i.e., accretion rate. We clearly see that to 
obtain the observed ratio the number density should be around 
2.6 - 8 x 10 4 cirr 3 . 

If the temperature is high enough, some methanol is re- 
leased in the gas phase. If «/ layers were formed on the grains 
and were all assumed to evaporate to the gas phase, the gas 
phase abundances of methanol would be given by 

n e iiisS , 
R = -5 = 3 x KTY-^ioWe. 

"ff 2 

Here, n g ~ n H2 \33 x 1CT 12 = «// 2 /-i2 the number density of 
grains, tih 2 is the number density of H2 in the gas phase, 5^0.3 
the fractional surface coverage of ;th species (e.g., methanol or 
water) in units of . 3 , S 6 the number of sites in the grain in units 
of 10 6 , and mo in units of 10 mono-layers. Since time taken for 
production of a mono-layer is ~ 0.1 Myr, by the time 10 mono- 
layers are produced, the grains are deeply inside in higher abun- 
dance region. Here the mono-layer production time is lower but 
the production efficiency of formation of methanol is lower as 
well. Thus R computed for methanol above, for «/ = 1 to 10 
varies from ~ 10~ 7 (s ifi . 3 ~ 0.03) to 3 x 10~ 5 (s ifi3 ~ 1), gen- 
erally agrees with the observational results. If, for instance, a 
fraction of grains are sublimated, say, /.12 ~ 0.1, the abun- 
dance could be even lower. For water, R ~ 10" 6 to 6 x 10~ 5 , 
where we have put «io = 1 to 10, Si$.3 = 0.33 for high abun- 
dance, and s ; o.3 = 2 for low abundance. This is also in the 
observed range. 

One of our important findings is that the parameter a(t) 
strongly depends on the population of the reactant species on 
the grain surface. This deviates significantly from unity. This 
seems to be a very important parameter, because in the usual 
rate equation we assume that a(t) is always 1 (a consequence 
of the assumption that the recombination is totally a random 
walk process). This is an overestimation. In Paper-I we also 
computed this parameter for the H2 molecule. We defined an- 
other parameter y called the 'catalytic capacity' and found that 
it goes down with the increased surface population. This shows 
that the rate of production indeed increases as the grain is fill- 
ing up. We also found that the behaviors of a(t) and y strongly 
depends upon the grain temperature. 

In our present calculation we restricted ourselves in two 
ways: (a) We considered only dense clouds where the accre- 
tion is such that most of the hydrogen is already in the molec- 
ular form in the gas phase. Thus accreting gas composition 
produces primarily water and methanol as is truly the case for 
dense clouds. In diffused clouds, on the other hand, the com- 
position is such that mostly molecular hydrogen is formed and 
are desorbed into the gas phase. Such a work was presented 



in Paper-I. (b) We restricted ourselves up to the formation of 
a mono-layer due to the fact that the complexity of the prob- 
lem rises with layer number and the unavailability of sufficient 
data (e.g., the binding energies at different surfaces). In reality, 
multi-layers would be produced and each layer is expected to 
have a different abundance due to the freezing out effect. For 
instance, in the first layers would be dominated by H2O and 
methanol, etc. In the later stage, accretion species will be dom- 
inated by CO, O2, N2, etc. Detailed results are in progress and 
will be reported elsewhere. 

The paper has been greatly improved due to the helpful 
comments of the anonymous referee who is acknowledged. The 
work of AD was supported by a RESPOND grant from ISRO. 
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Table 5. Surface coverage of the major species when one mono 
layer is built 



Major 


Used 


Ei, & Ed from Table 2 


Ej, & Ej from Table 3 


Ej, & Ed from Table 3 


species 


abundances 


Tunneling not allowed (Model 1 ) 


Tunneling not allowed (Model 2) 


Tunneling allowed (Model 3) 






CI 


C2 


C3 


C4 


CI 


C2 


C3 


C4 


CI 


C2 


C3 


C4 


O 


high 


4.98 


32.97 


33.29 





1.8 


13.13 


12.74 





1.72 


13.17 


12.83 





(in %) 


low 





12.56 


12.61 








13.82 


10.24 








13.67 


10.68 





o 2 


high 


19.32 








22.49 


25.26 


19.77 


19.76 


28.0 


25.39 


19.32 


19.43 


28.16 


(in %) 


low 


0.24 





(1 


0.38 


0.67 


0.39 


0.76 


1.05 


0.56 


0.35 


0.56 


0.6 


OH 


high 


8.92 


11.13 


10.97 





2.89 


2.35 


2.32 





2.92 


2.19 


2.34 





(in %) 


low 





0.2 


0.22 








0.15 


0.27 








0.15 


0.2 





H 2 


high 


16.33 


11.15 


11.19 


32.97 


13.36 


10.57 


11.23 


18.82 


13.23 


11.34 


11.17 


18.84 


(in %) 


low 


61.04 


48.59 


48.71 


61.43 


60.41 


46.81 


49.69 


60.27 


60.55 


47.02 


49.25 


60.47 


CO 


high 


50.45 


44.75 


44.54 


44.54 


48.60 


44.99 


45.38 


44.19 


48.61 


45.42 


45.69 


42.47 


(in %) 


low 


38.58 


38.57 


38.38 


38.1 


12.66 


12.41 


9.67 


12.31 


12.16 


12.68 


9.72 


12.12 


H 2 CO 


high 














3.86 


4.07 


3.86 


4.46 


3.96 


3.75 


4.07 


4.65 


(in %) 


low 


0.13 


0.07 


0.07 


0.09 


0.16 


0.15 


0.15 


0.15 


0.17 


0.12 


0.12 


0.16 


CH3OH 


high 














0.54 


0.64 


0.53 


1.06 


0.58 


0.53 


0.53 


1.16 


(in %) 


low 














26.01 


26.16 


29.03 


26.14 


26.41 


25.87 


27.2 


26.26 


C0 2 


high 














3.0 


1.29 


1.15 


3.47 


2.97 


1.3 


1.03 


4.72 


(in %) 


low 














0.08 


0.06 


0.11 


0.08 


0.14 


0.05 


2.18 


0.39 



Table 6. Time taken to build one mono-layer 





Used 
accretion 


Ej, & Ed from Table 2 
Tunneling not allowed (Model 1) 


E b & E d from Table 3 
Tunneling not allowed (Model 2) 


Ej, & Ed from Table 3 
Tunneling allowed (Model 3) 






CI 


C2 


C3 


C4 


CI 


C2 


C3 


C4 


CI 


C2 


C3 


C4 


Time 
in Year 


high 


1.6(3) 


1.3(3) 


1.3(4) 


3.6(4) 


1.7(3) 


1.6(3) 


8.3(3) 


3.2(4) 


1.7(3) 


1.6(3) 


5.7(3) 


3(4) 


low 


1.1(5) 


1.1(5) 


9.4(5) 


1.1(6) 


1.1(5) 


1.1(5) 


2.4(5) 


1.0(6) 


1.1(5) 


1.1(5) 


3.2(5) 


1(6) 



Table 7. Calculated value of a, for different species, using low 
abundances of the accreting species for Model 2 and keeping 
the grain at 10K. These simulations were carried out for the C2 
scheme by considering an Olivine grain having 10 4 number of 
sites. 







Surface coverage in (%) 


Species 


Site 


10 


20 


30 


40 


50 


60 


70 


80 


90 


100 


OH 


10" 


2.92 


3.08 


3.32 


3.53 


3.68 


3.79 


3.88 


3.98 


4.05 


4.12 


5 x 10" 


2.71 


2.78 


2.98 


3.12 


3.23 


3.32 


3.4 


3.47 


3.53 


3.59 


10 5 


2.61 


2.69 


2.83 


2.98 


3.09 


3.18 


3.25 


3.31 


3.37 


3.42 


H 2 


10" 


2.95 


2.97 


3.05 


3.27 


3.34 


3.46 


3.47 


3.48 


3.48 


3.47 


5 x 10" 


2.67 


2.69 


2.76 


2.92 


3.04 


3.1 


3.12 


3.13 


3.11 


3.1 


10 5 


2.58 


2.6 


2.68 


2.81 


2.92 


2.97 


2.994 


2.997 


2.997 


2.994 


o 2 


10 4 


3.24 


3.34 


3.61 


3.83 


4.03 


4.17 


4.3 


4.43 


4.55 


4.67 


5 x 10" 


2.55 


2.64 


3.04 


3.29 


3.5 


3.66 


3.8 


3.92 


4.03 


4.12 


10 5 


2.36 


2.52 


2.79 


3.07 


3.3 


3.46 


3.6 


3.71 


3.81 


3.89 


HCO 


10 4 


1.23 


1.24 


1.28 


1.47 


1.64 


1.78 


1.89 


2.04 


2.25 


2.36 


5 x 10" 


1.08 


1.14 


1.29 


1.43 


1.55 


1.68 


1.78 


1.89 


2.04 


2.18 


10 3 


1.07 


1.11 


1.25 


1.41 


1.54 


1.64 


1.74 


1.85 


1.98 


2.14 


H 2 CO 


10" 


3.25 


3.25 


3.25 


3.29 


3.4 


3.47 


3.49 


3.55 


3.71 


3.77 


5 x 10" 


2.77 


2.78 


2.81 


2.95 


3.06 


3.11 


3.13 


3.15 


3.22 


3.32 


10 5 


2.64 


2.66 


2.73 


2.84 


2.95 


3.0 


3.02 


3.06 


3.13 


3.24 


H 3 CO 


10" 


1.39 


1.38 


1.4 


1.44 


1.55 


1.64 


1.69 


1.81 


1.99 


2.06 


5 x 10" 


1.2 


1.21 


1.24 


1.37 


1.49 


1.57 


1.63 


1.69 


1.79 


1.89 


HP 


1.17 


1.19 


1.26 


1.37 


1.47 


1.54 


1.59 


1.65 


1.75 


1.88 


CH3OH 


10" 


3.33 


3.32 


3.33 


3.39 


3.49 


3.54 


3.56 


3.66 


3.79 


3.864 


5 x 10" 


2.86 


2.87 


2.9 


3.01 


3.12 


3.18 


3.21 


3.24 


3.33 


3.43 


HP 


2.71 


2.72 


2.78 


2.87 


3.0 


3.04 


3.08 


3.11 


3.19 


3.31 
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